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Abstract 

We reconsider quarkonium production in a field-theoretical setting and 
we show that the lowest-order mechanism for heavy-quarkonium produc- 
tion receives in general contributions from two different cuts. The first one 
corresponds to the usual colour-singlet mechanism. The second one has 
not been considered so far. We treat it in a gauge- invariant manner, and 
introduce new 4-point vertices, suggestive of the colour-octet mechanism. 
These new objects enable us to go beyond the static approximation. We 
show that the contribution of the new cut can be as large as the usual 
colour-singlet mechanism at high Pt for J ftp. In the ip' case, theoretical 
uncertainties are shown to be large and agreement with data is possible. 

Keywords: heavy-quarkonium production, vector-meson production, gauge in- 
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1 Introduction 



Years after the first disagreement between data JI!,'2J and tfie colour-singlet model 
(CSM) IHlll], the problem of heavy-quarkonium production -in particular at the 
Tevatron- is still with us. Indeed, it is widely accepted now that fragmentation 
processes, through the colour-octet mechanism (COM) 5J dominate the produc- 
tion of heavy quarkonia in high-energy hadronic collisions, even for Pt as low as 
6 GeV. In the COM, one then parametrises the non-perturbative transition from 
octet to singlet by unknown matrix elements, which are determined to repro- 
duce the data. However, fragmentation processes are known to produce mostly 
transversely-polarised vector mesons at large transverse momentum jHj and, in 
the J/ip and ip' cases, this seems in contradiction with the measurements from 
CDF [T. For a comprehensive review on the subject, the reader may refer to [H]. 

We therefore reconsider the basis of quarkonium (Q) production in field the- 
ory, and concentrate on J/ip and ip' production. We shall see that new contri- 
butions are present in the lowest-order diagrams, and we shall also explain how 
one can build a consistent and systematic scheme to go beyond the static ap- 
proximation. To this end, we shall use 3-point vertices depending on the relative 
momentum of the constituent quarks and normalised to the leptonic width of the 
meson. We shall show that, in order to preserve gauge invariance, it is required 
to introduce vertices more complicated than the 3-point vertex. 

Finally, we shall see that our formalism can be easily applied to the production 
of excited states. In the case of the theoretical uncertainties are unexpectedly 
large and allow agreement with the data. 

2 Bound states in QCD 

All the information needed to study processes involving bound states, such as de- 
cay and production mechanisms, can be parametrised by vertex functions, which 
describe the coupling of the bound state to its constituents and contain the infor- 
mation about its size, the amplitude of probability for given quark configurations 
and the normalisation of its wave function. In the case of heavy quarkonia, 
the situation simplifies as they can be approximated by their lowest Fock state, 
made of a heavy quark and an antiquark, combined to obtain the proper quantum 
numbers. Furthermore, it has been shown [Hj that, for light vector mesons, the 
dominant projection operator is and we expect this to hold even better for 
heavy vector mesons as this approximation gets better in the case of (j){ss). 

The transition qq ^ Q can then be described by the following 3-point func- 
tion: 

T^;\p,P) = Tip,Ph,, (1) 
with P = pi — P2 the total momentum of the bound state, and p = {pi + P2)/2 
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Figure 1: Bound-state vertex obtained by multiplying a point vertex, representing 
a structureless particle, by a form factor. 

the relative momentum of the bound quarks, as drawn in Figure H This choice 
amounts to describing the vector meson as a massive photon with a non-local 
coupling. 

We do not assume that the quarks are on-shell : their momentum distribution 
comes from r{p,P) and from their propagators. In order to make contact with 
wave functions, and also to simplify calculations, we assume that T{p, P) can be 
taken as a function of the square of the relative cm. 3-momentum p of the quarks, 
which can be written in a Lorentz invariant form as p ^ = —p^ + '"'^'j^} . For the 
functional form of r(p, P), we neglect possible cuts, and choose two otherwise 
extreme scenarios: a dipolar form which decreases slowly with p, and a Gaussian 
form: 



both with a normalisation N and a size parameter A, which can be obtained 
from relativistic quark models [TU]. We shall see in Section 4 how we fix the 
latter using the leptonic-decay width. 

3 Lowest-order production diagrams 

In high-energy hadronic collisions, quarkonia are produced at small x, where 
protons are mainly made of gluons. Hence gluon fusion is the main production 
mechanism. In the case of J/ip and if)', a final-state gluon emission is required 
to conserve C parity. This gluon also provides the Q with transverse momentum 
Pt- We assume here that we can use coUinear factorisation to describe the initial 
gluons, and hence that the final-state gluon emission is the unique source of Pt- 
All the relevant diagrams for the lowest-order gluon-initiated production process 
can be obtained from that of Figure |21 by crossing. There are six of them. 

These diagrams have discontinuities, which generate their imaginary parts. 
In order to find these, we can use the Landau equations [TTj. It is sufficient to 



r(p,P) 




and r(p,P) = Ne~^, 



(2) 
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Figure 2: Box diagram. 

consider the diagram of Figure El for which the Landau equations become: 

\i{cl-m^) = 

X2{{ci-Py-m') = 

Xsiic, - - m') = (3) 

A4((ci -h- - m^) = 

AiCi + A2(ci-P) + A3(ci-A;i) + A4(ci-A;i-A;2) = 

with m the quark mass. These equations have only two solutions in the physical 
region: one which is always present and gives a cut which starts at s = (A;i+fc2)^ = 
4m^, shown in Figure El (a), corresponding to a cut through the two s-channel 
propagators, and another one when the meson mass M is larger or equal to 
2m (see Figure |21 (b): this corresponds to a cut through the two propagators 
touching the meson). The latter leads to the colour- singlet model [3], which 
assumes that the quarks should be put on-shell to make the meson. The former 
cut has not been considered so far for the description of inclusive production. Let 
us mention however that similar cuts are dominant in diffractive production of 
vector mesons [121, DVCS [TBj . 

We are going to consider this s-channel cut in detail. To avoid complica- 
tions, we choose quark masses m > M/2 high enough for the second cut not to 
contribute. This will also simplify the normalisation procedure for the vertex. 
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(a) cut 1 



(b) cut 2 



Figure 3: The first family (a) has 4 diagrams and the second family (b) 6 diagrams 
contributing the discontinuity of gg — ^Sig at LO in QCD. 



3.1 Non- locality and gauge invariance 

The first problem one immediately faces when evaluating the diagrams of 
Figure El (a) is that of gauge invariance: whereas these diagrams are gauge in- 
variant if we have a photon instead of a Q, they are not for a finite-size object. 
Indeed, the vertex function r{p,P) takes different values in diagrams where ei- 
ther the on-shell quark or the antiquark touches the Q: the relative momentum 
p is then either p = 2ci — P or p = 2c2 + P, so that the delicate cancellation that 
ensures current conservation is spoilt. 




(a) (b) 



Figure 4: Illustration of the necessity of a 4-point vertex. 



The reason is easily understood: if one considers a local vertex, then the gluon 
can only couple to the quarks that enter the vertex. For a non-local vertex, it is 
possible for the gluon to connect to the quark or gluon lines inside the vertex, as 
shown in Figure |^. These contributions must generate a 4-point qqQg vertex, 

^In the cases of Figure |21 (b) , for M — 2ra, no gluon emission is kinematically allowed and 



5 



r^D(ci, C2, g, P). In general, its form is unknown, but it must obey some general 
constraints jEi [1^] : 

• it must restore gauge invariance: its addition to the amplitude must lead 
to current conservation at the gluon vertex; 

• it must obey crossing symmetry (or invariance by C conjugation) which 
can be written 

rifj(ci, C2, g, P, m) = -7orJfJ(-C2, -ci, q, P, -m)1'7o; (4) 

• it must not introduce new singularities absent from the propagators or from 
r(p, P), hence it can only have denominators proportional to (ci — P)^ — 
or (c2 + PY — m?- 



(3) 

• it must vanish in the case of a local vertex F/i oc 7^, hence we multiply it 

byr(2ci-P,P)-r(2c2 + P,P). 

These conditions are all fulfilled by the following choice JS] : 

rW(ci,C2,P,g) = -tgsT^A^i^c^ + P,P)-T{2c2-P,P)] 

7m (5) 



X 



(C2 + P)^ — (Ci — P)^ — TV? 

where the indices of the colour matrix T are defined in Figure and Qs is the 
strong coupling constant. 




Figure 5: The gauge-invariance restoring vertex, F*^^). 

It must be noted first that the quark pair (ci, C2) that makes the meson is now 
in a colour-octet state. We thus recover the necessity of such configurations, in 
this case to restore gauge invariance. We must also point out that this choice of 
vertex is not unique. We postpone a full study of the 4-point vertex ambiguity [T3] 
to another paper. 

When taken into account in the calculation of the discontinuity of gg —>■ Qg, 
T^^^ introduces two new diagrams^ shown in Figure IHl Including these contribu- 
tions in the calculation of the amplitude, we obtain a gauge-invariant quantity. 

the diagrams are directly gauge- invariant. 

^The contributions of the triple-gluon vertex on the left and r'"*^ on the right is zero due to 
charge conjugation. 
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Figure 6: New contributions to the discontinuity of the amphtude from T^'^\ 



4 Amplitudes 

We define A'^'^'"^ , i = 1,..., 6, as the unintegrated amphtude given by the usual 
Feynman rules for the four diagrams of Figure El (a) and the two of Figure IHl We 
choose the loop momentum i so that ci = £ + ki and C2 = £ — /c2- We then have 
for the imaginary part of the physical amplitude Ai: 

X 27r(5+((£ + kif - m2)27rr ((£ - fca)^ - m^) 

(6) 

This polarised partonic amplitude is thus obtained by contracting A'^'^'^" with 
the polarisation vectors of the gluons e\ j^, e\^^ and with that of the vector 
meson £3^^, by integrating on the internal phase space restricted by the cutting 
rules and by summing the six contributions from the diagrams of Figure El (a) 
and Figure ini 

The complete expressions of the polarisation vectors are as follows. One of the 
transverse polarisations can be taken as orthogonal to the plane of the collision: 

^1 —^2 —^3 —^4 — ^T, {() 

with eT-ki = eT-k2 = ex-P = 0. The other transverse polarisation can be taken 
as 



- ■ ^ i (^^1 + uk2 + sq) = el^ (8) 



for the two gluons, and as 



(10) 

for the final-state gluon and for the vector meson, where s = {ki + /c2)^, i = 
{k2 — qY and u = {ki — qY are the Mandelstam variables for the partonic process. 
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Finally, the longitudinal vector-meson polarisation can be taken as 

1 /. . fs + M^^ 



M 



ki + k2 



(11) 

(12) 



To complete the calculation of the amplitude, we need to normalise the 3- 
point vertex function F^^-*. We use here the leptonic decay width to fix this 
normalisation I16j. 

The width in terms of the decay amplitude is given by 



(13) 



where d2{PS) is the two-particle phase space [TTJ. 

The amplitude is obtained as usual through Feynman rules. At lowest order 
in as, only the 3-point vertex function needs to be considered (giving an explicitly 
gauge invariant answer). The square of the amplitude is then obtained from the 
diagram drawn in Figure [7| 




Figure 7: Feynman diagram for '^S\ — > l^. 
In terms of the sub-amplitudes A^''^ , B'^'^ and C^'^ defined in Figure [71 we have^: 
I \M\'d2{PS) = -l^.^'A^'^^^B^'^'^-^C^y, (14) 

where the factor A^^ = {g^^, — ^-j^) = Ylii^'i,iJ.^XiJ.' comes from the sum over 
polarisations of the meson and the factor 1 from the averaging over the initial 
polarisations. 

, after integration on the two-particle phase space, is found to be 



B 



pp 



3 



pppp 



M2 



(15) 



App' 



^We performed the calculation in the Feynman gauge, but the results are gauge invariant. 
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A^'^ (or equivalently C^'^'^) can be written (see Figure |H1): 



r(A;,P) 



(27r)4^ ^ ' ((A; - f )2 - m2 + ie)((A; + f )2 - + ie) ' 
with en the heavy-quark charge. 



(16) 




Figure 8: Feynman diagram for ^Si 



r 



Performing the integration by residues, one obtains 



with 



/(A, M, m) = 



-NI{A, M,m)A 



'^dKK'^V{k,P) {2K^ + ?,m^) 



(17) 



and K = \k\. J is a function of A through the vertex function T(k, P) and is not 
in general computable analytically, but it is straightforward to get its numerical 
value. 




1.4 1.6 1.8 

A (GeV) 



m.= 1 .60 GeV 
m^= 1 .69 GeV 



mc=1.78 GeV 
m^=1.87 GeV 
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1.4 1.6 1.8 

A (GeV) 
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m.= 1 .60 GeV 
m^= 1 .69 GeV 



=1.78 GeV 
=1.87 GeV 



Figure 9: Normalisation for a dipole (resp. Gaussian) form of T{p, P) in the J/ip 
case as a function of A: left (resp. right). 

We can then put all the pieces together using Eqs. (fT!?jl and (fT^ to determine 
N from the measured leptonic width. We show in Figure IHl the result in the J/ip 
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case. As can be seen, the normalisation of the vertex depends rather strongly 
on rric and A, but very little on the assumed functional dependence of the vertex 
function. We shall see later that once is determined from the leptonic rate, 
the production cross section depends little on these uncertainties. 



5 Production cross sections 

We can now evaluate the production cross section from the s-channel cut. As 
stated before, we assume that collinear factorisation can be used, in which case 
the link between the partonic and the hadronic cross sections is given by the 
following general formula: 

d^CT s d(7 

^T53= dxidx2 gixi) gix2) --^6{s + i + u- M^), (19) 
Jq tt at 

where xi and X2 are the momentum fractions of the incoming gluons, P = {E, P) 
is the momentum of the meson in the cm. frame of the colliding hadrons, g{x) 
is the gluon distribution function^ taken at the scale -y/ M ^ + P^ . 

In the cm. frame of the colliding hadrons, introducing the rapidity y = 
tanh~^(^) and the transverse momentum Pt, we obtain the double-differential 
cross section in Pt and y from Eq. (fTUI) : 

da /"^ , , , . , ^ ^„ da 



dydPj 



f da 

I dxidx2g{xi)g{x2)2sPT^S{s + i + u- M^) . (20) 
Jo dt 



At this stage, we can perform the integration on X2 (or xi) using the delta 
function. 

In terms of the transverse energy Et= ^P^ + M2, we get 

s = SX1X2, i = M"^ - xie"^v^ET, u = - xsc^Vs^t, (21) 
so that we obtain 

_ x,ErV-se-y - 

~ V~s{V~sxi - Erevy ^''^ 
The double differential cross section on P^ and y then takes the following form: 

da _ ^^ 2sPTg{xi)g{x2{xi)) da ^^^^ 



dydPr Jx^^^ \/s{\/sxi — Exey) di 
where x^*" corresponds to a;2 = 1 in Eq 



*In our calculations, we have chosen two LO gluon parametrisations, MRST jl8| and 
CTEQ j^. For each plot, the one used will be specified. 
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The last step is to relate the partonic differential cross section ^ to the 
amplitude calculated from our model. To this end, we use the well-known formula: 

daP'^''' 1 



,2 i« -^sl-^'^'^T, (25) 
dt IGiTs^ 

where |AlP'J''*p is the squared polarised partonic amplitude for gg Qg, averaged 
only over colour for polarised cross sections, and where p, q, r and s are the 
helicities of the four particles. 

As we are concerned with polarisation only for the Q, we sum over gluon 
polarisations and define, for r = L,Ti,T2: 

di ^ ^ di ' 

p,q,s=Ti,T2 

Finally , we have the double-differential polarised cross section on Pt and y : 
dar _ 2sPTg{xi)g{x2{x{)) dar ^27) 



dydPr J^min \^{y/sxi — Erey) di 



6 Results for J/ijj 

Before presenting our results, we need to choose a value for A and rric. Several 
studies have shown, in the context of relativistic quark models [10 , that the scale 
of the vertex function is between 1.42 and 2.6 GeV, and that is between 1.42 
and 1.87 GeV. We choose here a value of A in the middle range, A = 1.8 GeV 
(we shall see that small variations do not affect our results much) , and a value of 
rUc equals to the mass, mc= 1.87 GeV, in order to have a coherent treatment 
of all stable charmonium states. 

Setting ^/s to 1800 GeV and considering the cross section in the pseudorapid- 
ity range \ri\ < 0.6, we get the following results for J/ip production at CDF. The 
first plot (see Figure ITUl (a)) shows our result {cttot, o't and ai) for m = 1.87 
GeV and A = 1.8 GeV. 

These new contributions are compared with the usual LO GSM 
It must be stressed that (in the Feynman gauge) the main contribution comes 
from the 3-point function (Eq. (P)). As can be seen from Figure El (b), the 
term that restores gauge invariance (Eq. ^) contributes little: the square of its 
amplitude (see Figured (b)) is about 10 times smaller than the square of the 
amplitude containing only a 3-point vertex (see Figured (a)). Furthermore, the 
interference term between the diagrams of Figure 0] (a) and that of Figure E](b) is 
negative, so that the effect of Eq. © is to reduce the total amplitude squared. In 
Figures we show that the normalisation of the results using the decay width 
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has removed most dependence on the choice of parameters^. Interestingly, as 
figure 11 (b) shows, the dependence on A is neghgible once values of the order of 
1.4 GeV are taken. 




1e-08 [ — ' — ' — ' — ' — ' — ' — ' — ^ t — ' — ' — ' — ' — ' — ' — ' — ^ 
4 6 8 10 12 14 16 18 20 4 6 8 10 12 14 16 18 20 

Pt (GeV) Pj (GeV) 



(a) (Ttot, (Jt and ctl (b) aroT and crpert, 

Figure 10: (a) Polarised (ctt and ct/,) and total {(Ttot) cross sections obtained 
with a Gaussian vertex function, m = 1.87 GeV, A = 1.8 GeV and the MRST 
gluon distribution, to be compared with LO GSM. (b) The total (gauge invariant) 
contribution (plain curve) compared with that of the 3-point vertex of diagrams 
of Figure 111(a) (dashed curve) in the Feynman gauge (GTEQ). 

We see that the contribution of the new cut matters at large Pt- Noteworthily, 
it is much fiatter in Pt than the LO GSM, and its polarisation is mostly longi- 
tudinal. This could have been expected as scalar products of El with momenta 
in the loop will give an extra \/S contribution, or equivalently an extra Pt power 
in the amplitude, compared to scalar products involving Et- One can show that, 
for A = 1.8 GeV, = 1.87 GeV and for MRST structure functions, the longi- 
tudinal cross section falls as whereas the transverse cross section behaves 
asymptotically as l/p}p'^. The asymptotic power of pt changes by 5 % for A 
varying between 1.4 and 2.2 GeV. 

Recall that in our calculation the LO GSM is zero because M < 2m. For 
M > 2m, we should have added our contribution to that of the LO GSM at the 
amplitude level. The net result would then be fiatter and larger. However, it is 
clear that the enhancement factor would not be large enough to reach agreement 
with the data. 

^Instead of a factor 100 of difference expected from (^ Na-I I ^ have less than a factor 2 
at Ft = 4 GeV and a factor 3 at Pt = 20 GeV. 
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(a) Vertex function 



(b) A 



(c) A and rric 



Figure 11: (a) Comparison between the polarised cross sections obtained with 
the dipole and the Gaussian vertex functions (MRST); (b) Variation of the total 
cross section due to a change in A for a fixed value of the quark mass (MRST); 
(c) Variation of the total cross section due to a change in m and A (CTEQ). 



7 Results for ip' 

Although the normalisation to the leptonic width removes most of the ambiguities 
in the J/ip case, it is not so for radially excited states, such as the ip'. Indeed, in 
this case, the vertex function must have a node. We expect it to appear through 

a pre-factor, 1 multiplying the IS* vertex function. Explicitly, T2sip,P), 

for a node anode should be well parametrised by 



anode J (1 + ^)2 V anode J 



e^. (28) 

(1 + f^^^ ' " ' 

In order to determine the node position in momentum space, we can fix anode 
from its known value in position space, e.g. from potential studies, and take the 
Fourier transform of the wave function. However, the position of the node is not 
very well-known, and it is unclear how to relate our vertex with off-shell quarks 
to an on-shell non-relativistic wave function. The most remarkable thing is that, 
because the integrand has a zero, the integral / of Eq. (fT8|l entering the decay 
width calculation can vanish for a certain value anode = ao, which turns out to 
be close to the estimated value of the zero in the wave function. Because of 
their different momentum dependence, the integrals that control the production 
are not zero for anode = ao. Hence our normalisation procedure can in principle 
produce an infinite answer. Of course, this means that one cannot be at anode = oq 
exactly. However, if one is close to it, then it becomes possible to produce a large 
normalisation. Hence in the ip' case, our procedure can produce agreement with 
the data at low Pt- 
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Figure 1121 shows that for anode = 1-334 GeV, one obtains a good fit to CDF 
data at moderate Pt (note that the slopes are quite similar. This is at odds 
with what is commonly assumed since fragmentation processes -with a typical 
1/Pt behaviour- can also describe the data). The ip' is predicted to be mostly 
longitudinal. 



1 




e 1e-05 r 

T3 



1e-06 r : 

1e-07 ' ' ' ' ' ' ^ 

6 8 10 12 14 16 18 20 

Pt (GeV) 

Figure 12: Polarised [a^ and ai) and total {(Jtot) cross sections for ip' ob- 
tained with a Gaussian vertex functions, anode = 1-334 GeV, m = 1.87 GeV, 
A = 1.8 GeV and the CTEQ gluon distribution, to be compared with LO GSM, 
GSM fragmentation [4^ and the data from GDF pOj. 

This effect of the node in the ip' vertex function could also solve the p — tv 
puzzle as suggested in ^^Ij since a slight modification in the integrand of / can 
produce a large suppression in the p — n decay of the ip'- Such a modification 
in the integrand is indeed expected to come from the presence of an off-shell u 
in the p — tc decay instead of an off-shell photon for the leptonic decay. On the 
other hand, in the case of J/ip, no important effect are expected. 

8 Conclusion and outlook 

In this letter, we have shown that there are two singularities in the box dia- 
gram contribution to quarkonium production. We have chosen quark masses so 
that only the s-channel singularity (which is usually neglected) contributes, and 
vertices without cuts. 

On the theoretical side, we have begun to map the ingredients needed to go 
beyond the static approximation M = 2m. They involve the introduction of new 
4-point vertices that restore gauge invariance. We postpone a full study of these 
to a later pubhcation pUj . 
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On the phenomenological side, we have shown that in the J/ip case, this 
new singularity produces results comparable to those of the lowest-order GSM, 
i. e. too small to accommodate the Tevatron data. In the ip' case, ambiguities in 
the position of the node of the vertex function can lead to an enhancement, and 
to an agreement with the data. Hence it is not clear that the same mechanism 
has to be at work for IS* and 25" mesons. 

Our approach can be used in the case of P waves, where we would simply input 
a suitable 3-point vertex instead of taking higher derivatives of the wave function 
and of the perturbative amplitude, but also in a fc^-factorisation framework [21] 
(which would enhance the cross section in the J/ip case), and can be combined 
with contribution from COM fragmentation. Because the quarkonium are mostly 
longitudinal in this work and transverse in fragmentation, it seems possible to 
reach agreement with polarisation data. 
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